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This  paper  presents  a  high  efficiency  initial  rotor  position  estimation  method 
for  Permanent  Magnet  Synchronous  Motors  (PMSM).  The  approach  uses  the 
viable  inductance  model  to  analyze  the  optimal  motor  injection  sine  wave 
frequency  as  the  motor’s  test  signal.  Unlike  other  high-frequency  injection 
methods,  this  approach  does  not  require  trial  and  error  experiments.  The 
injection  frequency  is  identified  by  programmable  simulations  using  Matlab. 
Experimental  evaluation  of  3 -phase  PMSM  showed  that  this  injection 
frequency  optimization  method  works  successfully.  The  proposed  method 
can  find  the  optimal  injection  frequency  without  experimentation,  and 
outperforms  other  test  signals  in  terms  of  accuracy,  vibration  quantity  and 
noise.  This  microprocessor-developed  3 -phase  PMSM  control  driver  can  be 
applied  to  electrical  appliances,  machine  tools  and  automation. 
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1.  INTRODUCTION 

Sensorless  Permanent-Magnet  (PM)  motors  are  commonly  used  in  industrial  machinery  and  home 
appliances.  Their  advantages  include  a  simple  structure,  good  reliability,  high  power  density,  good  efficiency 
and  low  maintenance  costs.  In  addition,  such  motors  will  not  be  damaged  by  arcing  or  friction  at  high-speed 
rotation.  PM  motors  use  non-contact  phase  commutation  instead  of  mechanical  phase-changing  devices  such 
as  brushes.  To  activate  and  rotate  the  PM  motor  smoothly,  position  sensors  are  needed  to  detect  the  rotor 
position.  Position  sensors  can  be  photo-  or  magnet-based,  but  rotor  position  detection  may  be  costly  and 
unreliable. 

Implementing  motor  control  without  position  sensors  can  enhance  the  reliability  and  reduce  the  cost 
of  motor  sensors.  In  a  motor  control  system,  rotor  position  must  be  known  to  obtain  better  control 
performance.  In  recent  years,  several  methods  have  been  developed  to  estimate  rotor  position,  including  the 
counter  electromotive  method  [l]-[4],  the  extended  Kalman  filter  [5], [6],  the  artificial  neural  networks  [7], 
Finite  Element  method  [8], [9]  and  sliding  mode  current  observer  [10] -[12].  These  methods  are  based  on 
modern  control  theory  and  realized  through  the  use  of  microprocessors. 

However,  these  methods  can  only  work  when  the  rotor  is  rotating.  As  noted  by  Haque  [13],  when 
the  rotor  is  still,  measuring  the  initial  position  of  the  rotor  is  difficult  in  several  respects:  1)  A  still  motor 
produces  no  electrical-field  or  magnetic -field  change  which  can  be  measured  to  evaluate  rotor-position,  2) 
position  estimation  is  load  dependent  [14],  and  3)  it  entails  high  calculation  complexity  and  sophisticated 
hardware  requirements  [15]. 

In  1996,  Matsui  [16]  reviewed  five  methods  of  estimating  the  initial  position  of  sensorless  PM 
brushless  DC  motor,  including:  1)  Auxiliary  sensors,  2)  open  loop  controls,  3)  specific  gate  patterns,  4) 
arbitrary  starting  and  5)  salient-pole  motors.  Compared  with  Methods  2)  -  4),  Method  5)  starts  the  motor 


Journal  homepage:  http://iaesjournal.com/online/index.php/IJPEDS 


IJPEDS 


ISSN:  2088-8694 


n  157 


more  smoothly  and  prevents  reverse  rotation.  Method  5)  is  mainly  based  on  the  magnetic  saturation  of  the 
magnetic  conductive  material.  The  test  signal,  voltage  or  current  are  applied  to  the  motor,  and  then  the 
difference  of  the  corresponding  current  peak  value  is  measured  to  calculate  the  current  position  and  rotation 
speed  of  the  rotor  [17].  In  practice,  to  maintain  the  rotor  position  at  standstill,  the  energy  within  the  motor 
must  be  offset.  Once  the  test  signals  are  input,  an  identical  counter-directional  signal  must  be  input  to  zero 
out  the  net  torque.  This  signal  injection  method  is  referred  to  as  the  Salient  Pole  method. 

The  Salient  Pole  method  can  avoid  undesirable  events  such  as  rotor  back-kicking  or  rotor  oscillation 
prior  to  smooth  running,  making  this  method  superior  to  the  other  four  methods.  Based  on  the  Salient  Pole 
method,  several  approaches  are  provided  to  further  enhance  the  precision  of  initial-position  estimation. 
In  research  on  the  Salient  Pole  method,  Schmidt  et  al.  [18]  established  a  table  of  input- voltages  and 
corresponding  peak  currents  of  a  motor.  When  the  motor  starts,  a  three-phase  voltage  is  applied  to  it.  The 
induced  peak  current  is  then  measured  and  the  initial  position  of  the  rotor  can  be  obtained  based  on  the  table. 
Tursini  et  al.  [19]  found  that  evaluations  of  rotor  position  based  only  on  peak  current  are  subject  to  several 
uncertainties.  They  proposed  an  improved  method  which  applies  the  same  voltage  waveform  16  to  measure 
the  feedback  current,  using  fuzzy  theory  to  determine  rotor  position,  resulting  in  estimations  with  errors  of 
less  than  20  degrees. 

Hatakeyama  [20]  evaluated  the  rotor  position  by  applying  two  sine  waves.  The  amplitude  of  the  first 
high-frequency  sine  wave  is  small,  so  the  rotor  remains  still.  However,  the  amplitude  of  the  second  sine  wave 
is  high,  providing  sufficiently  strong  magnetic  saturation.  Since  the  first  sine  wave  cannot  induce  a  strong 
current,  the  detected  rotor  position  cannot  be  determined  to  be  in  the  one-to-two  quadrant  or  the  three-to-four 
quadrant.  Hence  the  second  sine  wave  is  applied  to  determine  the  quadrant  of  the  rotor  position.  However, 
the  strong  second  sine  wave  will  result  in  some  rotor  rotation. 


2.  RESEARCH  METHOD 

The  line  inductance  of  a  PMSM,  denoted  by  Lt ,  can  be  described  as  f  =Leq  +Lend  where  Leq  is  the  2- 
D  equivalent  inductance  and  Lend  is  the  end-turn  inductance  of  windings.  Lend  is  caused  by  three-dimensional 
(3-D)  magnetic  flux  in  a  PMSM,  and  is  usually  much  less  than  Lend  .  Thus,  this  paper  is  concerned  only  with 
Leq .  Some  assumptions  are  listed  here  to  simplify  the  analysis. 

The  line  inductance  of  a  PMSM,  denoted  by  Lt ,  can  be  described  as  f  =Leq  +Lend  where  Leq  is  the  2- 
D  equivalent  inductance  and  Lend  is  the  end-turn  inductance  of  windings.  Lend  is  caused  by  three-dimensional 
(3-D)  magnetic  flux  in  a  PMSM,  and  is  usually  much  less  than  Lend  .  Thus,  this  paper  is  concerned  only  with 
Leq .  Some  assumptions  are  listed  here  to  simplify  the  analysis. 

Al.  We  ignore  the  eddy  current  effect  on  inductance. 

A2.  The  shape  of  the  shoes,  the  teeth  and  the  slots  of  the  stator  is  simplified  as  rectangular  shapes. 

A3.  For  the  calculation  of  the  air-gap  permeance,  the  field  pattern  assumes  that  the  flux  lines  are  always 
orthogonal  to  the  intersection  surface  of  the  air-gap  and  the  steel. 

A4.  The  B-H  curve  of  the  PM  is  linear  with  the  constants  of  the  relative  permeability  ( jur )  and  the 
remanence  ( Br ). 

A5.  The  relative  permeability  of  the  rotor  and  the  stator  back  iron  are  infinite  (juFe=co)’  and  the  nonlinear 
interaction  of  these  two  fields  is  considered  only  to  calculate  the  air-gap  inductance. 

A6.  The  input  coil  current  is  constant. 

As  shown  in  Figure  1,  in  this  paper,  all  three  phases  are  excited  by  different  winding  currents,  as 
opposed  to  the  approach  used  in  [21].  It  is  known  that  Leq  is  the  sum  of  the  air-gap  (Lg)  and  the  slot 

inductances  ( Lslot ). 

4 q=(Lg+Lslot)^  (1) 

B  where  Ns  is  the  total  number  of  slots,  and  the  denominator  3  is  the  number  of  the  phases  of  a 
PMSM.  The  inductances  Lg  and  Lslot  are  respectively  functions  of  the  flux  linkage  and  . 
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Figure  1.  Simplified  shoes,  teeth,  and  slots  structure 


Figure  2  shows  the  flux  linkages  of  the  windings  divided  into  five  parts:  , 

where  is  the  flux  crossing  the  air-gap,  while  consists  of  ,  ®s2  ,  and  . 


5  and 


Figure  2.  Flux  pattern  of  winding  current 


The  Lslot  can  also  be  divided  into  four  parts  as: 

Lslot  =  24l  +  Ls2  +  24-3  +  Ls4 


(2) 


Figure  3.  (a)  equivalent  magnetic  circuit  of  ®5l  .  (b)  equivalent  magnetic  circuit  of  <&s2  .  (c) 
Equivalent  magnetic  circuit  of  air-gap  flux 


Figure  3(a)  and  (b)  respectively  show  the  equivalent  magnet  circuits  of  ®5land®52  •  Note  that  N  is 
the  number  of  turns  of  each  stator  tooth,  Rsl  is  the  magnetic  reluctance  of  ®5l ,  and  Rs2  is  the  magnetic 
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reluctance  of  052  .  It  follows  from  Figure  3  that  05l  =  3M  /  2 Rsl .  According  to  the  definition  of  the  inductance, 
we  obtain 


3ATJOd  9N 


,vl  . 

di  4  R, 


\sl 


(3) 


It  is  well  known  that Rsl  1  ^]^)Z  ,  where  and  are  defined  in  Figure  1  The  current  direction  of 
faces  b  and  c  are  the  same,  therefore  ®s2  and  Ls2  are  almost  equal  to  zero: 

^2  ~  9  Bs2  ~  0 

In  Figure  2,  i.e.,  the  slot  between  the  teeth  of  phases  b  and  c.  The  magnetic  motive  force  (MMF)  at 
the  top  of  the  slot  is  equal  to  2Nt ,  while  it  is  zero  at  the  bottom.  This  is  because  the  turns  are  uniformly 
distributed  in  the  slot,  and  the  magnetic  field  Hs 3(x)  crossing  the  slot  increases  linearly  from  the  bottom  of 
the  slot  to  the  top.  Since  Hs3(0)  =  0  and  Hs3(h)  =  2Ni/ (b0  +  T-t) ,  we  have  Hs3(x)  =  3Ni/2[(b()  +  T-t)h\ ,  where  *  is 
the  distance  from  the  bottom  of  the  slot  (see  Figure  2). 

This  allows  us  to  express  43  from  the  energy  viewpoint  as 


3  A)jV2fe 

^  i2  4  (fio  +  T-t) 


(5) 


It  is  well  known  that  ^''4  ~  Ls2 ,  therefore  and  Ls4 


are  almost  equal  to  zero: 


0^4  « 0  Ls 4  «  0  ^ 

The  equivalent  magnetic  circuit  in  Figure  3(c)  describes  the  air-gap  flux  ®gthat  goes  from  the  shoe 
of  phase  a  to  the  reluctance  of  phase  a  RFa  and  the  PM  reluctance  Rg  ,  then  1/2  <4>g  separately  pass  through 
reluctance  of  phase  a  and  phase  b.  RFa  ,  RFb  ,  RFc  and  Rg  can  be  approximated  as: 


_  _|_  § 
MrMOTZ  MOTZ 


RFa  = 
RFb  = 
RFc  = 


fr)  |  h 

MrMaABa\)^  MoMal^al)* 

to  !  h 

HrHb\(Bb\)TZ  /Jo/Ub2(Bb2)tZ 

hp  ,  h 

Mr  Mel  ( Bc\  ycz  MoMcl  ( Bc2  )tz 


(7) 


Mai ,  Mb l  ?  Mc i  are  relative  degrees  of  permeability  in  the  tooth,  and  ^a2  ,  ^bl ,  ^c2  are  those  in  the 

shoe.  It  then  follows  from  Figure  3(c)  that  Og  =  3M  /  2[Rg  +  RFa  +  Rss  ] 

^s:s  =(^  +RFb)(Rg  +RFc)  I  (Rg  +RFb)  +  (Rg  +RFc )  •  T^US 


_3¥^°g_9  A2 
2  di  4  [/?£  +  +  7?sv  ] 


(8) 


However,  by  assumption  A5,  the  relative  permeability  of  the  stator  back  iron  is  a  nonlinear  function 
of  the  flux  densities  in  the  shoes  and  the  teeth  of  the  stator.  These  flux  densities  are  defined  as  follows: 
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(9) 


where  Om  are  the  fluxes  of  the  PM  flowing  through  the  three  phase.  are  the  equivalent  slot  linkage  fluxes 
of  the  three  phases. 

It  is  easy  to  calculate  ,  ®mb  ,  ,  ®sa  ®sb  ®sc  with  the  following  equations: 


®ma  Z,\  Bmgadl  ®mb  Z,\  Bmgb^  ®mc  ^  J  Bmgc^ 


0  =hL  0  : 

^sa  2  ]\f 


AN 


l 

0~ 


4i 

4Af 


(10) 


It  is  apparent  that  flux  densities  are  functions  of  ,  so  the  flux  0g  should  make  the  following 
error  function  £(0g)  zero  from  the  equivalent  circuit  shown  in  Figure  3  (c): 


£(<D?)  =  |M-<Dg[ffg  +RFa  +i?J»0 

Z  (11) 

£(®g)  is  a  nonlinear  equation  with  a  single  variable  ®g  .  Applying  the  iteration  method  to  the  above 
equation,  we  can  easily  obtain  0g  that  makes  £(®g)«0,  provided  that  ®m  is  known. 

We  substitute  the  flux  0^  into  (8)  to  calculate  the  air-gap  equivalent  inductance  Lg .  The  slot 
inductance  Lslot  and  the  air  gap  inductance  Lg  can  be  substituted  into  Eq.  1  to  calculate  the  magnetic 
saturation  and  the  motor’s  effective  equivalent  inductance  Leq .  Finally,  substitute  into  (1)  to  obtain  Leq 

Figure  1  shows  the  d-q  winding  model  of  a  PMSM.  In  this  figure,  Id  and  lq  represent  the  armature 
d-q  windings,  2d  the  constant  field  winding  corresponding  to  the  permanent  magnet  mounted  on  the  rotor  and 
3d  and  3q  the  equivalent  windings  for  eddy  current.  In  this  paper,  the  iron  loss  arises  from  the  resistance  of 
the  equivalent  winding  for  eddy  current  [1]. 


3.  ROTOR  POSITION  ESTIMATION  METHOD 

As  noted  in  [2],  the  motor  inductance  changes  with  the  rotor  angle,  thus  the  relationship  between  the 
inductance  value  and  the  current  can  be  written  as: 


When  the  voltage  is  kept  constant  in  the  above  formula,  the  inductance  value  is  inversely 
proportional  to  the  current  differential,  so  the  rotor  position  can  be  determined  by  calculating  the  difference 
of  the  inductance  value  from  the  speed  of  the  current  response. 

As  shown  in  Figure  4,  the  drive  circuit  of  the  three-phase  permanent  magnet  synchronous  motor  is  a 
three  half-bridge  circuit.  We  use  a  three-phase  leg  with  a  total  of  six  transistors  and  flywheel  diodes  (  S{  and 
S2  ,  S3  and  S4  ,  S5  and  S6 )  to  form  a  three-phase  inverter,  where  va  is  the  input  voltage  of  phase  a,  vb  is  the 
input  voltage  of  phase  b,  vc  is  the  input  voltage  of  phase  c,  the  neutral  point  voltage  at  the  intersection  of  the 
three  phases  is  vn  and  Vdc  is  the  DC  voltage  source.  The  three-phase  permanent  magnet  synchronous  motor 
coil  Y  passes  phase  a  as  phase  current  ia  ,  phase  b  as  phase  current  ib  ,  and  phase  c  as  phase  current  ic  . 
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Figure  4.  Three-phase  motor  and  drive  circuit 


The  three-phase  voltage  input  method  includes  six  modes,  where  the  a-phase  injection  definite 
voltage  signal  mode  is  taken  as  an  example: 


vab  ^an  ^bn)  vac  (—  van  v cn )  Vdc 

ia  =  ia(+)  k  =  h  =  “0 1 %(+) 

% 

where  Vdc  is  the  signal  amplitude 

vab  =Vdc  =  2  Rsia(+)  + 1  (L  + 

Similarly,  assuming  the  a-phase  injection  is  a  negative  voltage  signal  mode 


(13) 


(14) 


The  positive  value’s  square  wave  test  signal  Vdc  is  input  to  the  motor  model,  the  simplified 

a 

resistance  is  =/?  and  the  signal  is  the  positive  inductance  value. 

Thus  we  can  rewrite  (13)  as  a  voltage  signal  relationship 

vab  =  Vdc  =  ia(+)R  +  Reqa(+)  ^  ^  ^ 


d-Vr/ra«)  ,ra(+)  =4 qa(+)IR 

The  negative  value’s  square  wave  test  signal  -Vdc  is  input  into  the  motor  model.  The  negative  signal 
input  uses  the  negative  inductance  value  ;  similarly,  the  voltage  is 


vab  Vdc  ia(-)R  +  Reqa(- ) 


dia(-) 

dt 


(18) 


The  current  response  equation  is: 

ia(-)(0  =  ^a-e~t,Ta(~h  ,ra(_}  =4 


(19) 
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We  set  a  current  sampling  time  tp ,  the  negative  current  \-)(tp)  and  positive  current  . 

Subtracting  the  absolute  values  from  each,  the  three-phase  current  difference  A i  is: 


A4=l4(-)(^)|-|4(+)(^)l 

A 4  ^  4(-)  )  I —  I  ic(+)  (t^ )  I  (20) 

Through  this  method,  we  can  use  the  magnitude  relationship  of  the  three  A i  values  to  separate  the 
360  degree  electrical  angle  into  six  parts,  each  with  an  electrical  angle  resolution  of  sixty  degrees,  as  shown 
in  Figure  5  In  an  electrical  angle  of  30  degrees  A ib  and  A ic  rise  and  fall  nearly  linearly,  making  this  area  a 
dividing  point,  and  the  other  five  blocks  are  treated  with  the  same  logic. 


Figure  5.  Electrical  angle  pardoning 


This  square-wave  method  can  be  used  to  successfully  identify  the  initial  rotor  position.  However,  if 
the  motor  has  a  large  stator  coil  resistance  value,  the  current  response  will  be  little  changed,  resulting  in  an 
excessively  large  judgment  error  in  determining  the  initial  rotor  position.  Therefore,  this  paper  seeks  to  select 
an  appropriate  injection  frequency  to  maximize  the  current  distance  and  minimize  the  rotor  position 
estimation  error. 


4.  OPTIMUM  FREQUENCY  SEARCH 

The  proposed  optimum  frequency  search  method  is  described  in  following  flow  chart.  Figure  6 
illustrates  the  process  of  determining  the  optimal  injection  frequency.  The  first  step  is  to  input  the  motor’s 
geometric  parameters.  Step  2  sets  the  coil’s  initial  current  value.  Of  particular  note,  the  coil  current  i{n)  is  not 

arbitrary,  but  is  rather  determined  iteratively.  The  current  iteration  calculation  seeks  to  determine  the  optimal 
current  ia(f optimal)  as  we  approach  the  optimal  frequency  foptimai ,  and  thus  calculate  the  circuit  current  i{n)  at 

the  time  of  analysis  to  optimize  the  accuracy  of  current  value  calculations. 

Step  3  inputs  the  data  produced  in  Steps  1  and  2  to  the  model  derived  in  Section  2  to  calculate  the 
motor’s  equivalent  inductance. 

Step  4:  The  optimal  injection  frequency  search  method  inputs  different  frequency  sine  wave  signals 
(from  1Hz  to  500Hz).  When  the  voltage  is  positive,  L^qa^(0e)  is  the  inductance  value;  when  y^  is 

negative,  4^(_)(<%)  is  the  inductance  value,  and  this  signal  is  input  into  the  Simulink  motor  model. 
Subtracting  the  absolute  value  of  the  half-cycle  coil  current  y+)  from  the  negative  half  cycle  coil  current 
obtains  A i ,  and  the  frequency  at  which  the  maximum  current  max(A/a)  occurs  is  the  optimal  current  foptimai 

Finally,  compare  the  optimal  frequency  current  and  the  current  at  the  time.  If  they  are  not  similar, 
then  add  a  current  error  ierror  and  continue  with  further  iterations.  If  they  are  similar,  terminate  the  operation 
with  the  optimal  frequency  simulation.  Where  the  current  error  value  ierror  is  (iopt  -\n))l  2,  if  they  are  almost 
equal,  the  operation  terminates  to  complete  the  optimal  frequency  search. 
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Figure  6.  Flow  chart  of  optimum  frequency  search 


Figure  7  shows  the  result  of  the  frequency  search  from  1Hz  to  500Hz.  The  largest  A i  occurs  at  21Hz, 
followed  by  20.05mA. 


400  SOD 


Figure  7.  Result  of  frequency  search 


5.  EXPERIMENTAL  RESULTS 

Experiments  used  three-phase,  16-pole,  and  9-slot  PMSMs,  with  microcontrollers  to  produce  the 
injection  signal  and  measure  the  coil  current  to  verify  the  simulation  results.  The  injection  signal  is  a 
sinusoidal  function  with  a  peak  value  of  75  V.  To  find  the  optimal  input  sign  wave,  the  frequency  of  the  input 
signal  was  increased  from  10  Hz  to  50  Hz.  The  relation  between  the  three-phase  current  difference  ( A i )  and 
the  rotation  angle  is  shown  in  Figure  8  The  amplitude  of  the  three-phase  current  difference  ( A i )  is  maximized 
when  the  frequency  of  input  sin  wave  is  20  Hz.  This  result  is  identical  to  the  simulation  result. 
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Figure  8.  The  relation  between  feedback  current  difference  ( A i )  and  rotation  angle.  The  input  signal  is  a 
sinusoidal  function  with  a  frequency  of  (a)  10  Hz,  (b)  20  Hz,  (c)  30  Hz,  (d)  40  Hz,  and  (e)  50  Hz 


As  shown  in  Figure  5  rotor  partition  estimate  accuracy  increases  with  the  three-phase  current  ( A i ) 
error.  To  improve  estimate  resolution,  this  paper  uses  the  arctangent  function  to  reverse  the  initial  electrical 
angle.  This  method  uses  sum-to-product  formula.  We  already  know  that  Aia,Aib,Aic  each  have  120  degrees 
difference  of  electrical  angle.  Using  Eq.  21,  we  can  obtain  Aiu,Aiv,Aiw  and  eliminate  the  DC  offset. 


A iu  =  (A ia  -Aib)/  *]3 

<  Aiv=(Aib-Aic)/^3  (21) 

Aiw=(Aic-Aia)/s/3 


We  again  use  the  plot  difference  to  obtain  the  x,y  for  each  90  degree  difference  of  electrical  angle: 


x  =  Aiu 

y  =  (A/y  —  A/w)  /  y/3 


(22) 


Finally,  the  electrical  arctangent  function  is  used  to  obtain  the  initial  angle: 


X 

6e  =  arctan(— ) 

y 


(23) 


Figure  9  shows  the  relationship  between  the  actual  position  and  the  estimated  rotor  position  to 
evaluate  the  estimated  accuracy.  We  use  two  injection  frequencies  of  the  developed  method  to  compare  with 
the  two  sine  waves  method.  The  result  proves  that  our  Sensorless  Initial  Rotor  Position  Estimation  Method 
(Figure  9(a))  get  the  better  result  than  Two  Sine  Waves  Method  [20]  (Figure  9(c))  at  the  error  of  actual  rotor 
position  and  the  estimated  value.  The  experimental  RMS  error  is  also  calculated  in  Table  1,  which  tells  the 
frequency  of  20Hz  in  injection  can  get  the  most  accurate  rotor  position. 
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Figure  9.  Rotor  position  estimates  for  (a)  novel  method  with  20Hz  injection  frequence  (b)  novel  method  with 
50Hz  injection  frequence  (c)  two  sine  waves  method  [20] 


Table  1.  Comparison  of  Different  Injection  Frequence  and  Methods 


Novel  method  with  20Hz 

Novel  method  with  50Hz 

Two  sine  waves  method 

RMS  Error 

11.44 

18.95 

18.37 

6.  CONCLUSION 

In  this  paper,  we  analyze  the  methodology  of  Magnetic  Circuit  Analysis  to  establish  an  equivalent 
inductance  model  of  permanent-magnet  synchronous  motors.  With  this  model,  different  rotor  positions  are 
found  to  correspond  to  different  inductor  values.  Therefore,  The  optimum  injection  frequency  is  obtained  by 
the  simulation  process  of  the  Matlab.  The  RMS  error  of  the  experimental  shows  the  optimum  injection 
frequency (20Hz),  which  is  much  lower  than  the  other  injection  frequency  and  injection  method.  Our  method 
can  help  the  driver  developer  in  their  early  motor  developmrnt  process  to  avoid  using  trial  and  error  to  find 
the  optimum  injection  frequency.  Besides,  using  this  injection  method,  there  is  only  one  frequency  injection 
applied.  This  can  exempt  the  excessive  starting  of  the  motor  that  would  usually  generate  the  noise  and 
vibration.  We  conclude  that  this  novel  invention  can  be  widely  used  in  the  intelligent  appliances,  that  the 
permanent  magnet  synchronous  motor  is  one  of  the  most  important  components  in  the  device. 
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